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Abstract. This is an expository paper on Lyapunov stability of 
equilibria of autonomous Hamiltonian systems. Our aim is to clar- 
ify the concept of weak instability, namely instability without non- 
constant motions which have the equilibrium as limit point as time 
goes to minus infinity. This is done by means of some examples. 
In particular, we show that a weakly unstable equilibrium point 
can be stable for the linearized vector field. 

I. Introduction 

Stability of the equilibrium is a mathematical field more then two 
centuries old. Indeed, Lagrange stated the celebrated Lagrange-Dirich- 
let theorem in the eighteenth century, and some so called converses of 
that statement are still proved nowadays. So many mathematicians 
have been interested in stability that we refrain from mentioning them 
with the exception of the most important, Lyapunov, who defended 
his doctoral thesis "The general problem of the stability of motion" in 
1892. The applications are also countless in mechanics and in most 
sciences. To start with the rich literature on this matter, see Arnold et 
al. pp, Meyer et al. [4], and Rouche et al. [7]. 

Important mathematical objects related to the instability of the equi- 
librium are asymptotic motions. Before their formal definition, let us 
mention that the upper position of a simple pendulum, and zero veloc- 
ity, constitute an unstable equilibrium and its asymptotic motions are 
neither rotations (when the pendulum swings around and around) nor 
librations (when it swings back and forth), and they stay between the 
two behaviors. 

Let us consider a smooth vector field / on an open set A C R N with 
an equilibrium point x G A, so f(x) = 0. We say that : (—00, b) — > A 
is an asymptotic motion in the past to the equilibrium point x, if 4>(t) 
is a non-constant solution to the o.d.e. x = f{x) such that (f)(t) — > x as 
t — > —00. In the sequel we briefly write 'asymptotic motion' instead of 
'asymptotic motion in the past' since we are only concerned with this 
kind of asymptotic motions. Of course the existence of an asymptotic 
motion implies the Lyapunov instability of the equilibrium point. The 
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basic sufficient condition for the existence of an asymptotic motion is 
the presence of an eigenvalue of f'(x) with strictly positive real part, 
see for instance Hartman [3] remark to Corollary 6.1, p. 243. 

In this paper we focus on autonomous Hamiltonian systems so in the 
sequel N = 2n, x = (q,p), q,p G M n , and the vector field is 

(1) (d p H(q,p),-d q H(q,p)) 

for some smooth H called the Hamiltonian function. Our aim is to 
clarify the concept of Lyapunov instability without asymptotic motions 
that we briefly call weak instability. This is done by means of some 
examples. 

Section [2] deals with linear systems. Of course there is a trivial situa- 
tion where weak instability appears: the free particle. The equilibrium 
is non-isolated and the eigenvalues vanish, the example can be done in 
one degree of freedom so in dimension 2. A more subtle instability of 
the equilibrium for a linear system is obtained when the eigenvalues are 
purely imaginary and some Jordan blocks have dimension greater than 
one, of course this can happen only in dimension at least 4. The exam- 
ple we are going to see comes from the planar restricted 3-body problem 
at one of the relative equilibria, the Lagrange equilateral points, also 
called the Trojan points, at the critical Routh value of the mass ratio 
of the primaries. 

In Section [3] we move on nonlinear systems. Their equilibria can 
be unstable even if we have stability for the linearized system as the 
Cherry Hamiltonian in dimension 4 shows by means of an asymptotic 
motion. Cherry's system is the third example of this paper, it was 
published in 1925 and, in the last 20 years, it became important in 
plasma physics, see Pfirsch [6] and the references therein. 

Our fourth example, also in dimension 4, comes from [TT] and shows 
that we can have weak instability of an Hamiltonian equilibrium which 
is linearly stable. Some systems, produced by Barone-Netto and my- 
self [9] and [TO] , preceded [UJ, they give non-Hamiltonian examples of 
weak instability for linearly stable equilibria. 

Hopefully, the concept of weak instability will stimulate further re- 
searches in stability within mathematical physics, together with other 
fresh notions like the "weak asymptotic stability" introduced by Or- 
tega, Planas-Bielsa and Ratiu, see [5] and the references therein. 



2. Weak instability for linear systems 

2.1. Free particle. Our first example is a particle on a straight line 
under no forces 
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(2) H (q,p) = ^, q,peR. 
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The Hamiltonian vector field is 

(3) (d p H(q,p),-d q H(q,p)) = (p,0). 

It is a linear field with the double eigenvalue 0. The integral curves are 

(4) q(t) = q(p)+p(0)t, p(t)=p(0). 

Each (g , 0) G IR 2 is an equilibrium point and its instability can be 
shown by means of the sequence (g(0),p(0)) = (q ,l/m) — » (qo,0) as 
m — > +oo. There are no asymptotic motions. 

2.2. Linearization at L 4 . Our second example is the quadratic part 
of the Hamiltonian function of the planar restricted 3-body problem 
at one of the relative equilibria, the Lagrange libration point L4 at the 
critical Routh value of the mass ratio of the primaries. In the sequel 
Q = (Qi,Q2),P = (pi,P2), (q,p) = (91, 92, Pi, P2), and 

111 1 

(5) H(q,p) = -y=det(p,q) + -\q\ 2 = -y=(pi92 -P2?i) + +? 2 )- 

The Hamiltonian vector field is 

(d Pl H(q,p),d P2 H(q,p), -d qi H(q,p), -d q2 H(q,p)) = 

(6) 



(q 2 / V2, -qi/ V2, -qi + p%j V2, -q 2 - Pi 



0. 



see H and the o.d.e. at the end of p. 256, with £ = q, rj = p, u = 1/ \/2, 
5 = 1, and also H at p. 258 in Meyer et al. jl]. 

It is a linear vector field with the double eigenvalues A = ±z/\/2 

and Jordan blocks ^ ^\ . The origin is now the unique equilibrium 

point. 

The function \q\ 2 is a first integral. Suppose the integral curve 
(q(t),p(t)) — > as t — > —00, then |g(t)| 2 = and this fact further 
implies that |p(t)| 2 = 0, indeed for q(t) = we have 

(7) j t W)\ 2 = 2p(t) ■ ( K -q 1 (t)+p 2 (t)/V2,-q 2 (t)-p 1 (t)/V^j 

So the integral curve is constant and we do not have asymptotic motions 
to the equilibrium point. 

The origin is an unstable equilibrium point as we can see with 

gi(t) = ^ cos 71' g2(t) = -^ sin 71' 
w t t t t 

Pl (t) = --cos- 7 =, P2 (t) = -sin-, 

for w (gi(0),g 2 (0),p 1 (0),p 2 (0)) = (1/m, 0,0,0) -> as m +00. 

Incidentally, in connection with the nonlinear 3-body problem which 
has the Hamiltonian vector field defined by ^ as linearization at L4, 
the book jl] at the end of Sec. 13.6 says that in 1977 two papers claimed 
to have proved the stability of the equilibrium, however one proof is 
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Figure 1. Asymptotic motion for Cherry Hamiltonian 

wrong and the other is unconvincing. The last sentence is: "It would 
be interesting to give a correct proof of stability in this case, because 
the linearized system is not simple, and so the linearized equations are 
unstable" . 

3. Instability for linearly stable equilibria 

3.1. Cherry Hamiltonian. Next, the famous Cherry Hamiltonian 
system shows that the equilibrium can be unstable even if it is sta- 
ble for the linearized system, briefly even if it is linearly stable. In 
Cherry [2] p. 199, or in Whittaker [8] p. 412, we can see the Hamilton- 
ian function H : 1R 4 -> K 

(9) H{q, p) = X - {q\ + pf) - (q 2 2 + p 2 2 ) + a [q 2 {q\ - p 2 ) - 2q lPl p 2 S ) . 

The Hamiltonian vector field, written as a column vector, is 

( px - 2oq 2 p x - 2oq x p 2 \ 

-q%- 2aq 2 q 1 + 2ap x p 2 
\ 2q 2 + op\ - aq\ J 

The linearized vector field (pi, — 2p 2 , — qi, 2q 2 ) is obtained for a = 0. 
The origin is stable for the linearized systems which consists of two 
harmonic oscillators: q\ = —q%, q 2 = —^q 2 . The eigenvalues are dis- 
tinct ±z, ±2i. However, the origin is Lyapunov unstable for the vector 



field (10) whenever a^O since it has the following asymptotic motion 



defined for t < 



sint , . sin(2£) 

qi(t) = — =—, q 2 (t)- 



[11) 



V2at' " w 2at 
, . cost . . cos(2t) 
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3.2. Variation-like Hamiltonian. Our final example shows that the 
origin is an unstable equilibrium point which is linearly stable and has 
no asymptotic motions for the system defined by 

(12) H(q,p) = Pip 2 + qiq 2 + a q\q 2} a^O. 

It is a particular case of the following Hamiltonian function introduced 
in [TO 

(13) H(q,p) = p lP2 + g( qi )q 2 , g(0) = 0, g'(0) > 0, 

where g G C 1 on a neighborhood of 0. The Hamiltonian vector field is 

/ P2 \ / P2 \ 



(14) 



V 



P2 
Pi 

-g'(qi)q2 
-gin) ) 



o(\(q,p)\)- 



P2 
Pi 

-9'(0)q 2 

\-g'(o)qiJ 

The origin is stable for the linearized system which consists of two 
harmonic oscillators: q\ = —g'(0)qi, q 2 = —g'(0)q 2 . In this case the 
eigenvalues are double ±i^g'(0) however the Jordan blocks are one- 
dimensional. 

The subsystem of the first and last canonical equations 



(15) 



Qi = P2, 



P2 



-g(qi 



separates. If we take a solution (qi(t) , p 2 (t)) of this subsystem and 
plug qi(t) into the second and third canonical equations, we then get 
the equations of variation of (15) along the solution (qi(t) , p 2 {t)) . This 
is why the function in formula (13) is called variation-like Hamiltonian 
in the title of this subsection. 

There are no asymptotic motions, indeed if the solution 

(16) (qi(t),q 2 (t), Pl (t),p 2 (t)) ^0 as t ->■ 



OG 



then (qi(t),p 2 (t)) = 0, since the origin is a local center for (15), and 
this implies (q 2 (t) , pi{t)) = too. 

In spite of this fact, the origin is unstable for (14) for most functions 
g as above. Theorem 3.3 in [TT] proves that stability is equivalent 
to the isochrony of the periodic solutions of the subsystem (15) in a 
neighborhood of 6 K 2 , and this implies the isochronous periodicity 
of all integral curves of (14) in a neighborhood of G M 4 . Moreover, 
Corollary 2.3 in [TT] for a smooth g provides 

5s"(0) 2 



(17) 



<T(o) 



3<7'(0) 



as the simplest necessary condition for (local isochrony and then) sta- 
bility. So the choice g(q\) = qi + crqf of the Hamiltonian (12) gives 
instability for all a ^ 0. In Figure [2] we can see the projection on the 
qi j 92-plane of the integral curve of the Hamiltonian vector field given 
by ph. 
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Figure 2. Projection of an unbounded orbit for H = 



Finally, let us remark that the Hamiltonian (12), composed with the 
symplectic transformation (Q, P) \-t (Q 1 + Q 2 , Qi — Q2, Pi + P2, Pi — 
P 2 )/y/2, becomes 

(18) I {Qi + Pi) -\{Ql + Pi) + ^=(Qi + Q2) {Ql - Qt) 

a function with some features in common with Cherry's Hamilton- 
ian @. 
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